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Abstract 

> 

. We demonstrate numerically the existence of nonlinear Tamm oscillations at the interface be- 

\ tween a substrate and one-dimensional waveguide array with both cubic and saturable, self-focusing 

o 
o 



> 



and self-defocusing nonlinearity. Light is trapped in the vicinity of the boundary of the array due 
to the interplay between the repulsive edge potential and Bragg reflection inside the array. In 
the special case when this potential is linear these oscillations reduce themselves to surface Bloch 
oscillations. 
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The problem of surface waves which could exist at the interface between two different 
media has been studied for several decades. These waves are, for example, investigated at 
a metal- vacuum interface air- water interface at the interface between two dielectrics 
, and in finite anti-dot lattices |^]. Tamm was the first to consider neglected edge effects in 
a semi-infinite Kronig-Penney model P]. He has discovered that strongly localized surface 
states could appear providing that the surface potential perturbation is strong enough. 
However, due to various experimental difficulties, his discovery was awaiting more than half 
a century on the experimental confirmation in semiconductor superlattices [^. From the 
practical point of view, surface waves represent convenient tools to investigate the properties 
of material interfaces. Furthermore, they could be implemented in various devices such as 
bio-sensors [7] and polariton lasers p. 

Recently it has been suggested that these surface states may also exist at the interface be- 
tween a homogenous medium (substrate) and a nonlinear waveguide array (NWA) |^. NWA 
represent arrays of evanescently coupled optical waveguides. Up to date, NWA have been 

_|l2L quadratic Q 
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^ saturable 

Q, Q, breathers Q, 



Mm 

ICl, and Bloch oscillations 



successively fabricated in materials exhibiting cubic [1 
Q, and nonlocal nonlinearity Discrete solitons 

diffraction management ^|, modulational instability 

[igj l are just a few examples of phenomena which have been observed in such systems. 

Very recently, the first experimental observation of discrete surface solitons in AlGaAs 
NWA exhibiting a cubic self-focusing nonlinearity [2^ has triggered further investigations 
of surface waves at the interface between a NWA and a substrate. The existence of sur- 
face gap solitons in the lattice with cubic self-defocusing nonlinearity has been reported in 



2l|. A crossover from nonlinear surface states to discrete solitons was studied, too 



lef. 

2^ . In these two papers it has been revealed that the vicinity of the edge enables a stable 
pro pag ation of various localized modes, such as fiat-top modes and inter-site ones (mode 
B) j2J|. Very recently, stronglv localized surface waves have been independently observed 
in NWA exhibiting saturable [2^ |3| and quadratic nonlinearity 0] , respectively. The fact 
that a surface can support the existence of localized states has been, for example, further 
exploited in vector j23| , nonlocal j3l as well as binary NWA 2^ , and in a soliton array case 



In this Letter we reveal the existence of Tamm oscillations at the edge of a semi-infinite 
NWA. These oscillations are the result of an interplay between a repulsive potential which 



2 



originates from the boundary of the array and the array's periodicity. We calculate this 
potential for a few different media and reveal that Tamm oscillations are more likely to 
occur in systems with stronger coupling. In a special case when the repulsive potential is a 
linear function of the distance from the edge of the NWA Tamm oscillations reduce to the 



well-known Bloch oscillations 
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switching, amplification, and steering 



IS 



3c I , which have considerable potential in all-optical 



32|. 



The propagation of light in periodically stratified structures such as NWA can be de- 
scribed analytically either by the Floquet-Bloch model 3, or by a coupled-wave theory 



13L Il7| . The model equation within the first approach reads: 
dy 2k dz^ rig 



(1) 



The propagation coordinate is along the y-axis, the amplitude of the electrical field is denoted 
by while k = 27ms/ X represents the wave number. Here, A is the wavelength of the used 
light in vacuum, and Ug is the extraordinary refractive index of the substrate. The periodical 
modulation of the refractive index which defines the nonlinear WA is denoted by n{z) while 
Auni is the nonlinear refractive index change [Auni << ns). In media with cubic nonlinearity 
Arini = AuqI while in media with saturable nonlinearity we have An„; = AuqI /{I + Id), 
where / is the peak light intensity and Id is the so-called dark irradiance Q|. 

On the other hand, light propagation through NWA may be described, within the tight- 
binding approximation, by the following set of conveniently normalized, linearly coupled, 
nonlinear, ordinary differential equations: 

■dUn . ^,rr . rr r r N /Irr |2 



+ C{Un+l + Un-1 - 2Un) ~ g{\UnV)Un = 0, 



(2) 



where q{\Up\'^) = /3|f/nP in cubic (c) media and g{\Un\'^) = /3|f/nP/(l + l^nP) in saturable 



(s) media 



35l |. while the nonlinearity coefficient (3 < for the self- focusing (/) and (3 > 



for the self-defocusing (d) case, respectively. Here C is the coupling constant, ^ is the 
)ropagation coordinate and f/„ is the normalized electric field envelope in the n-th waveguide 
3^. Integral of motions are power and Hamiltonian: 

n 

H, = J2{C\Un-l-Un\^ - (3[\n{l + \Un\^) - \Un\^]}, 



n-l 



(3) 



Assuming stationary solutions of staggered form [/„ = F„exp[2(— z/^ + nn)] {y represents 
soliton frequency) for defocusing cases, and of unstaggered form f/„ = F„exp(— zi/^) for 
focusing cases, together with the assumption |-Fo| ^ \F±i\ ^ I -^±2 1 for on-site (A) mode and 
|-^±i| ^ l-^±2| ^ \F±?\ for inter-site (B) mode, one may find the following expressions for 
the maximal amplitude in the array: 



In addition for A mode -Fn(>o) = 01^ Fq and for B mode -Fn(>i) = tt"'^^-^!, in both saturable 
and cubic media. Here a = -^C/{v — 2C) for (d) and (/) case, respectively. Examples of 
the oscillatory behavior of both on-site mode (mode A) and mode B are depicted in Fig. 1. 
Here the energy is too low to overcome the repulsive potential from the edge of the array 
so both modes start to move away from the interface until they are back-reflected because 
of the Bragg condition. As this traversing is usually accompanied with radiation reflected 
modes loose power and eventually do not reach back to the flrst channel of the array. As a 
result, consequent oscillations have longer and longer periods and modes gradually run away 
from the edge. The trapping of localized modes at one channel after one or more oscillations 
is possible as well. Here is interesting to mention that surface waves could exist due to 
the balance between self-bending and deflection from the edge of a bulk photorefractive 



crystal, too [37]. Equivalently, the behavior of moving localized modes can be interpreted 
by the interplay between a repulsive force (due to the boundary effects) and the effective 
Peierls-Nabarro potential (due to system discreteness) j^^. 

We calculate the effective repulsive potential as the difference between the truncated 
potential and a reference potential with periodic boundary condition, and present results 
for symmetric localized structures centered on-site in an array that consists of 13 elements. 
As the peak of this localized structure (with the highest amplitude Fq) approaches the edge 
more and more terms of the corresponding Hamiltonian have to be truncated. For example, 
if the peak of mode A resides in the 6th channel one can obtain the following expressions 
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FIG. 1: Tamm oscillations of a) mode A with C = 0.5, |/3| = 3.34, and soliton frequency u = 2.3 
launched into the first channel of the array, and b) mode B with C = 0.5, \f3\ = 3.34, and soliton 
frequency u = 4.83 launched into the second and third channel of the array. 

for the repulsive potential: 

V;W = C|Fo|V°-C|Fo|2(l + a)V°-0.5/?|Fo|V^ 



V}^1 = C|Fo|V°-C|Fop(l + a)V° 

-/3|Fo|V2 + /3/n(l + ai2|Fo|2), (7) 



in both saturable and cubic media in defocusing and focusing cases. After a straightforward 
calculation it is possible to get explicit expressions for Vrep, Vrep, etc. where the number in 
the superscript denotes the number of the truncated channels. For example, for the Tamm 
state in the self-focusing cubic case we have found: 

-C |Fop(l + a)2(aio + ^8 ^ ^6 ^ ^4 ^ ^2 ^ 
-0.5 /3\Fo\\a^^ + + a'^ + a'^ + + (8) 

Note that here (3 < and a = —C/{v — 2C). The dependence of the repulsive potential 
on the distance from the edge of the array is presented in Fig. 2. Stronger coupling (i.e. 
shorter coupling length) results in stronger repulsion (see Fig. 2a) while stronger nonlinearity 
decreases the repulsive potential (Fig. 2b). Beams which are strongly pushed off from the 
edge will experience Bragg reflection inside the array earlier than weakly rejected beams. 
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FIG. 2: Dependence of the effective repulsive potential on the distance from the interface between 
a substrate and the array, a) Self-focusing case: \f3\ = 3.1, \Fq\'^ = 0.24 and different values of C. 
b) Self-defocusing case: C = 0.5, |i*bP = 0.172, and different values of /3. 



Thus, the same input beam will have shorter spatial periods of Tamm oscillations in arrays 
with stronger coupling and weaker nonlinearity. 

The value of the corresponding repulsive potentials for both cubic and saturable, self- 
focusing and self-defocusing cases also depends on the soliton frequency u. For a fixed value 
of this frequency and fixed values of parameters C and (3, the following relation is usually 
fulfilled: V^^p > V^':^^ > V^^^ > Vf'/p but interior elements may permute their position as 
well. This relation may explain why this effect has not been observed in recent studies in 
AlGaAs waveguide arrays exhibiting a cubic self- focusing nonlinearity 0, . Interestingly, 
in the case when is a linear function from the interface distance the input beam will 
experience Bloch oscillations [19||, which have a considerable potential for application in 
all-optical devices as reported, for example, in Ref. 32. Restricting ourselves here only 
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FIG. 3: (Color online) Tamm oscillations in a waveguide array where channels are 4 //m wide and 
separated by 4.4 ;um while Ano = 3.7x 10~^. The input pattern with amplitude ratios of 1 : 0.5 : 0.1 
has been launched into the first channel, a) Self-defocusing saturable case An„; = 6.2 x 10^'*, 
r = I / III = 6, b) self-defocusing cubic case An„i = 4.42 x 10^'*, c) self-focusing saturable case 
An„i = 3.34 X 10~^, r = 6, and d) self-focusing cubic case An„i = 2.62 x 10~^. 

to NWA, this phe_n wa. studied in curved NWA M in NWA «t,t a hneariy 

growing effective refractive index [19]. Please note that, from the surface's point of view, 

'or li ght intensities lower 
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Tamm oscillations do not belong to nonlinear effects (i.e. they exist 
than the threshold for the onset of highly nonlinear Tamm states 2 > 

In order to check our findings we performed additional simulations based on a numerical 
solution of Eq. (1). We, arbitrarily, have used the parameters which are achievable in 
lithium niobate waveguide arrays exhibiting a self-defocusing saturable nonlinearity (Ang = 
3.7x 10-^ Auni < 0.001) and green light with A = 532 nm. For such samples the periodically 
modulated refractive index can be well approximated by a cos^ function Numerical 
results which have been obtained by virtue of a beam propagation method are given in 
Fig. 3. Note that in reality this effect occurs in a highly elongated space (20mmx0.1mm 
in Fig. 3a). 

As the position of the input beam shifts towards the interior of the array, the period 
of Tamm oscillations increases. Fig. 4 may be understood as a numerical proof that these 



oscillations are indeed a linear surface effect. Namely, the light intensity necessary to 



Tamm oscillations (Fig. 4a) is not high enough to form a surface soliton p, |2 
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FIG. 4: (Color online) Light propagation in a saturable nonlinear NWA with the same input 
profile and lattice parameters as in Fig. 3. Here Arini = 6.21 x 10~^ and r = 6. The input beam 
is launched into a) the first channel and b) into the third channel of a semi-infinite lattice. 



l25l . l26l | but suffices to form a narrow breather in a channel which is only two channels away 
from the substrate-array interface (Fig. 4b). 

Experimentally, Tamm oscillations may be observed, for example, in one-dimensional 
ithium niobate waveguide arrays, where sample lengths up to 50 mm have been achieved 



2i 



26l |. However, up to now we failed in a direct observation of Tamm oscillations from 



the top of the sample because of the rather low average level of scattered light that is 
superimposed by larger scattering amplitudes from small surface defects. An alternative 
approach has been outlined in Ref. [19|, and corresponding experiments may be performed 
in the near future. 

In conclusion, we have demonstrated numerically the existence of Tamm oscillations at 
the interface between a substrate and a one-dimensional homogeneous nonlinear waveguide 
array. Light is trapped in the vicinity of the edge of the array due to the interplay between the 
edge repulsion and Bragg reflection. Approximate analytical expressions for the repulsive 
truncated potential are given for different types of nonlinear interactions as well as for 
different system parameters. These oscillations reduce to the case of Bloch oscillations when 
the repulsive potential is a linearly decreasing function of the distance from the edge of the 
semi-infinite nonlinear waveguide array. 
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